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A system consisting of a spring, dashpot, and mass upon which is 
mounted an eccentric driven by a motor with a linear torque-speed char- 
acteristic, is analyzed by perturbation procedures based on small reciprocal 
of rotational inertia. Periodic solutions of the third order system, which 
arises when the angular position of eccentric mass is taken as the new 
independent variable, are constructed, and their stability is analyzed. An 
asymptotic solution is also obtained which is more general than a periodic 
solution, in that the averaged rotational speed is a slowly varying function, 
rather than a constant. The results are applicable to the determination of 
the interaction between the rotational motion of a flexibly mounted motor 
and the translational vibratory motion of its frame. 

I. INTRODUCTION 

In a recent paper Senator analyzed a system consisting of a spring, 
dashpot, and mass upon which is mounted a rotating eccentric weight 
driven by a motor with a linear torque-speed characteristic. 1 This 
system has been analyzed by several authors, under different as- 
sumptions on the values of the parameters of the system (see Refs. 5 
through 9), and Senator discusses their results. The system is a model 
for the interaction between the rotational motion of a motor driving 
an eccentric and the translational vibratory motion of the frame, which 
is caused by this rotation. 

In Ref. 1, Senator constructed periodic (rotational) solutions by 
means of a perturbation technique, based on small reciprocal of ro- 
tational inertia. However, he did not analyze the stability of the periodic 
solutions directly, but proceeded in a somewhat different manner. Thus, 
he introduced a van der Pol type transformation, but imposed a sub- 
sidiary condition on the slowly varying functions of time which differs 
from the one usually imposed in the method of averaging. He then 
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made assumptions regarding the order of smallness of various deriv- 
atives, and dropped all second order terms from the equations for the 
slowly varying quantities, obtaining what he called averaged equations. 
The stationary solutions of these averaged equations correspond to 
periodic solutions of the original system, and he analyzed the stability 
of the stationary solutions on the basis of the corresponding linearized 
variational equations. 

It is the purpose of this paper to show how the stability condition 
obtained by Senator may be derived rigorously for sufficiently large 
values of rotational inertia. This is done by taking the angular position 
of eccentric mass as the new independent variable, constructing periodic 
solutions of the resulting third order system, and then analyzing the 
linearized variational equations corresponding to them. Perturbation 
procedures, based on the small reciprocal of rotational inertia, are used. 

An asymptotic solution is also obtained which is more general than 
a periodic solution, in that the averaged rotational speed is a slowly 
varying function, rather than constant. However, this asymptotic 
solution is not completely general, in that the transients in the trans- 
lational motion, which decay on a much faster scale, are not included. 
The asymptotic solution nevertheless provides insight into the manner 
in which a stable periodic solution is approached, and the analytical 
results are borne out by some numerical calculations. 

II. PERIODIC SOLUTIONS 

The equations of motion, in dimensionless form, for the system under 
consideration are (from Ref. 1), 



d 2 u . nu du 



ddY . „ d 2 d 



— I sin 6 — -r-s cos 
dr/ dr 



(1) 



d 2 6 / , dd d 2 u \ , Q x 

" = Jp - b- a -7-2 cos ■ (2) 

dr Y dr dr I 

Here r is dimensionless time, a, b, p and f > are constants, and 
e > 0, the reciprocal of dimensionless inertia, is a small parameter. 
Also, u is the dimensionless translational displacement, and is the 
angular position of eccentric mass. Instead of dealing with the fourth 
order system (1) and (2), as did Senator, it turns out to be more con- 
venient to take as the new independent variable. Accordingly, denning 

= f (3) 

dr 
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the third order system 

Je i 1 + "* cos e fe) + erf cos 9 1P s=€( p- h ® ® 

is obtained. 

A periodic solution to (4) and (5) is sought in the form 

u = u(6, e) m u u (6) + eU,(0) + e 2 u 2 (d) + ■■■ (6) 

12 = £2(0, e) - co + <£2,(0) + e 2 £2 2 (0) + • • • (7) 

where w,(0) and £2,(0) are periodic in 0, with period 27r, and o> is a 
constant. Substitution of (6) and (7) into (4) and (5), and comparison 
of the lowest powers of e, leads to 

^ + 2ra^ + Wo = a^sin0 (8) 



Uii 



W ° ~tf0~ aa) ° "dP = ( p _ ° w o) • ( 9 ) 

The periodic solution of (8) is 

u = aw 2 A [(l - a»o) sin — 2fa> cos 0] (10) 

where 

a = [(i - <4) 2 + 4r 2 ^]- 1 . (id 

In order that £2j(0) should be periodic, it is necessary from (9) that 

p = bu + aw 2 / COS -~i) = bu + a'fcoo A = p*(u ) (12) 

using (10). ( ) av denotes average over a period 2tt of 0. Equation (12) 
gives a relationship between w„ and p, the dimensionless stall torque, 
and this relationship is depicted graphically in the figure for a = 0.707, 
f = 0.2, 6 = 0. It is noted that w is a triple valued function of p in 
part of the range. Senator concluded from his analysis that the middle 
branch corresponds to unstable periodic solutions, while the outer 
branches correspond to stable ones, a result verified in this paper. 
Now, from (9), (10), and (12) it follows that 

£2 t = {«i - }a s «gAo[(l - «j|) cos 20 + 2fco sin 20]} (13) 

where w, is a constant, which is to be determined from the condition 
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that fl 2 (0) should be periodic. It is clear as to how the higher order 
terms in the expansions in (6) and (7) may be obtained, but they will 
not be needed in the subsequent analysis. The periodic solutions in (6) 
and (7) are equivalent to those derived by Senator as periodic solutions 
of (1) and (2). It is necessary, of course, to perform a quadrature of 
equation (3) in order to obtain a relationship between and t. 

III. STABILITY ANALYSIS 

The variational equations corresponding to the periodic solution 
U, fi, given by (6) and (7), are formed by substituting 

t* = (w + £), Q = (fi + v) (14) 

in (4) and (5), and linearizing in £ and 77. Thus, 

^ + -g, + 2f (a| + f,) +f 
(1+— ft("S+8')+- fl 8—i 

+ 2e«fl cos 6 jp V + e«(2 2 cos 6 + ebr? = 0. (16) 

Equations (15) and (16) are linear equations with periodic coefficients, 
and the form of solution is known from Floquet theory. 2 Moreover, 
if all the characteristic exponents of the variational equations have 
negative real parts, then the periodic solution U, tt is asymptotically 
stable. The behavior of the characteristic exponents will be analyzed 

for < e « 1. 

The limiting case e — > 0+ will first be considered. In this case, from 
(6) and (7), Q = w and u = u {6) so that, from (15) and (16), d-q/dd = 
and 

* w + 2 ** % + { - 4»° 8i » e - - $> - f f^ (i7) 

Hence one of the characteristic exponents is X = 0, and the remaining 
two characteristic exponents satisfy 

(co Xo) 2 + 2f(co X ) + 1 = (18) 

and hence have negative real parts, since f > and « > 0. For sum- 
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ciently small e, these real parts will remain negative, so that it suffices 
to investigate the characteristic exponent which vanishes as e — » 0. 

In the light of Floquet theory, a solution of (15) and (16) is sought 
in the form 

£ = e xe P(d), n = e xe Q(6) (19) 

where P and Q are periodic in 0, with period 2r, and 

X = eX, + e 2 X 2 + • • • (20) 

P(6) = P (6) + eP,(0) + e 2 P 2 (6) + • • • (21) 

Q(e) = Qo(e) + eQ^e) + t 2 Q 2 (e) + • • ■ . (22) 

It is a straightforward matter to substitute from (6), (7), and (19)-(22) 
into (15) and (16), and to compare like powers of e. In particular, it 
is found from (16) that dQ /dd = 0. Omitting a multiplicative constant 
and taking Q = 1, it is then found that 



2 d2p ° _i_ o d 2 u dP Q du Q . 



and 

+ 2ao. cos 6 — § + c«4 COS ^ + 6 = 0. (24) 

Now, in order for Qi(6) to be periodic, it is necessary from (24) that 

cd^ + b + acoJ cos ^py = (25) 

where 

#o = (co P + 2u ) (26) 

is periodic in 0, with period 2t. But, from (23), 

co 2 f| a + 2fco S + Rq = 2 (aco 2 sin + f Wo f || + „ ) (27) 

and w is given by (10) and (11). Straightforward calculations lead to 
R Q = 2«oS A 2 {[(1 - co 2 ) 2 (2 - co 2 ) + 4«(1 - 2co 2 )] sin 

- fw [(l - co 2 )(5 - co 2 ) + 12^2] cos 0}. (28) 
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Thus, from (25), 

WoX, + b + a 2 & A 2 [(l - co 2 )(5 - co 2 ) + 12f 2 co 2 ] = 0. (29) 

However, as may be verified from (12), using (11), equation (29) 
may be written in the form 

coo*: + *£ = 0. (30) 

Since the sign of X in (20) is determined by the sign of \ u for suffi- 
ciently small e > 0, it follows that the periodic solution u, 9, is asymp- 
totically stable if dp*/da > 0, and is unstable if dp*/du> < 0. That is, 
the middle branch of the figure corresponds to unstable periodic solu- 
tions, while the outer branches correspond to asymptotically stable 
ones, provided that e > is sufficiently small. 

IV. MORE GENERAL SOLUTIONS 

In this section a more general solution of (4) and (5) is constructed, 
for < c <3C 1. Thus an asymptotic solution is sought in the form 

u = v (o>, d) + ^.(co, 6) + e 2 v 2 (co, 6) + • • • (31) 

ft = co + ew^u, 6) + e 2 u> 2 (co, 6) + • • • (32) 
where w,-(w, 0) and w,(co, 0) are periodic in 0, with period 2ir, and 

f = e 9l (co) + *%(*) + • • • • (33) 

This procedure may be regarded as a variant of the method of aver- 
aging. 3 The above solution is more general than the periodic solutions 
constructed previously, since the latter correspond to the case in which 
co is constant, rather than a slowly varying function of 0. However, 
this solution is not completely general, since the initial transients in 
(4) are not taken into account. 

Substituting (31) and (32) into (4) and (5), using (33), and comparing 
the lowest powers of e, it follows that 



co 2 y$ + 2fco fj + .o = «co 2 sin 6 (34) 

do ou 



and 



"{ft + g ) + ao}2 cos 9 w = (p ~ M - (35) 

The periodic solution of (34) is 
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v n = aco 2 A(co)[(l - w 2 ) sin 6 - 2fw cos 0] (36) 



where 



A(«) = [(1 - or) 2 + 4fV]- x . (37) 

In order that u\ should be periodic, it is necessary from (35), using 
(36), that 

cog i(co) = [p — bu — a 2 fo>°A(a>)]. (38) 

Then ir 1 (o), 0) may be found from (35), to within an arbitrary function 
of w. This arbitrariness is usual in averaging procedures, and may be 
removed by requiring (u\-(w, 0)} nv = 0, so that, from (32), w is the aver- 
aged value of Q. Higher order terms in the asymptotic expansion (31) 
and (32) may be obtained in a systematic manner. 
Now, from (33) and (38), 

^ = 6 [p - p»] + 0(0 (39) 

(lu a) 

where 

p*( W ) = 6w + aVA(co). (40) 

As previously remarked, the case in which co is constant corresponds 
to the periodic solutions constructed earlier. From (11), (12), (37), 
and (40), it follows that co is the lowest order approximation to a 
stationary solution of (39). If w ^ w , equation (39) determines, to 
lowest order, the slow variation of ta with 0. The direction in which co 
changes is determined, to lowest order in e, by the sign of [p — p*(«)], 
and is illustrated in Fig. 1 for p = 0.45, to which there correspond 
three values of co , denoted by w ol , o m , and w 0r • 

Under more general initial conditions similar results should hold, 
for sufficiently small e, provided that the initial value of £2 is not too 
close to coom • This is because Q does not change significantly, for suffi- 
ciently small e, during the time in which the initial transients in the 
translational motion die out. 

A partial check of these analytical results was made by Senator, 4 who 
carried out some numerical solutions of (1) and (2). With a = 0.707, 
f = 0.2, b = and e = 0.1, he chose initial conditions consistent with 
the unstable periodic solution corresponding to p = 0.425, that is, the 
periodic solution corresponding to p*(w 0m ) = 0.425. He then carried 
out numerical solutions of (1) and (2) for p = 0.45 and p = 0.4. He 
found that for p = 0.45 the solution approaches the periodic solution 
corresponding to co 0r in the figure, that is, to p*(co 0r ) = 0.45, while for 
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Fig. 1 — Stall torque vs. averaged rotational speed. 

p = 0.4 the solution approaches the periodic solution corresponding 
to p*(o3 i) = 0.4. These results are consistent with our analytical re- 
sults. Moreover, the number of cycles required before the solution 
settles down to the appropriate periodic solution was somewhat larger 
in the case p = 0.45 than in the case p = 0.4. This is consistent with 
(39), and the presence of the factor (1/w) multiplying [p — p*(co)] 
therein. 
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